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Abstract 

The recessive opaque-2 mutant gene (o2) reduces oc-zeins accumulation in maize endosperm, changes the amino acid 
composition of maize kernels, induces an opaque endosperm, and increases the lysine content of kernels. The quality 
protein maize (QPM) inbred line CA339 (o2o2) and an elite normal inbred line liao2345 (0202) were used to construct o2 
near-isogenic lines (NILs) by marker-assisted selection (MAS) using the co-dominant SSR marker phi057. Two specific o2 NILs 
were constructed, named liao2345/o2-1 and liao2345/o2-2. However, the kernel phenotypes of the two o2 NILs were 
different from each other. liao2345/o2-1 had the wild-type vitreous endosperm, which is similar to its recurrent parent 
liao2345, while the endosperm of liao2345/o2-2 was opaque, identical to typical o2 mutant individuals. In comparison to 
their recurrent parent liao2345, the lysine concentration of liao2345/o2-1 was similar and the lysine concentration in 
liao2345/o2-2 was doubled. SDS-PAGE analysis indicated that liao2345/o2-1 had the same zeins ratio as liao2345, whereas 
the zeins concentration of liao2345/o2-2 was markedly lower. Sequence and transcript abundance analyses indicated that 
the CDS of two o2 NILs are derived from CA339, but they have different promoters. The 02 transcript of liao2345/o2-2 is 
largely inhibited because of an rbg transposable element inserted between the TATA box and initiator codon of liao2345/ 
o2-2. We concluded that different crossing-over patterns during the process of o2 NIL construction resulted in the different 
kernel phenotypes of the two o2 NILs. We surmise that the reversion of liao2345/o2-1 to wild type was due to the 
recombination with the wild type liao2345 promoter during introgression and backcrossing. The o2 mutant gene of donor 
(CA339) is a null mutant because of low 02 expression. However, its CDS probably encodes a protein with normal function 
which can maintain the normal accumulation of zeins in maize endosperm. 

Citation: Chen Y, Zhou Z, Zhao G, Li X, Song L, et al. (2014) Transposable Element rbg Induces the Differential Expression of opaque-2 Mutant Gene in Two Maize 
o2 NILs Derived from the Same Inbred Line. PLoS ONE 9(1): e85159. doi:10.1371/journal.pone.0085159 

Editor: Tai Wang, Institute of Botany, Chinese Academy of Sciences, China 

Received July 6, 2013; Accepted November 24, 2013; Published January 9, 2014 

Copyright: © 2014 Chen et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits 
unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. 

Funding: This work was supported by the China Agriculture Research System (CARS-02-01), the Beijing Municipal Natural Science Foundation (6112003). The 
funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. 

Competing Interests: The authors have declared that no competing interests exist. 

* E-mail: limingshun@caas.cn (ML); zhangshihuang2013@163.com (SZ) 

9 These authors contributed equally to this work. 



Introduction 

Maize (£ea mays L.) is one of the most widely cultivated cereal 
crops in the world. Dependence on maize as a protein source puts 
people at risk of dietary protein deficiency because maize protein is 
deficient in two essential amino acids: lysine and tryptophan. 
Therefore, maize is a poor source of protein for humans and other 
monogastric animals [1]. In order to increase the quality of maize, 
a long-term selection experiment was initiated in 1896 by CG 
Hopkins at the University of Illinois [2] that resulted in the Illinois 
High Protein Strain (IHP). IHP has 2.5 times the amount of 
protein contained in normal maize and the majority of this 
increase is due to zeins, which contain no lysine [3]. In 1964, 
Mertz reported that the o2 mutant gene can alter the amino acid 
composition in maize endosperm and double the lysine content 
[4]- 

Because the o2 mutant gene can increase the lysine content in 
maize kernels, breeders began using this mutant to cultivate high- 



quality maize with more lysine. However this mutant maize had 
poor agronomic traits. This gene was cloned by transposon 
tagging [5,6] and its sequence was published [7]. In the 1990s, the 
mechanism through which the 02 gene regulates protein 
expression was intensively studied. The 02 gene encodes a protein 
that has structural homologies to transcriptional activators [8]. 
The OPAQUE2 protein has a "leucine-zipper" motif that can 
bind to zein DNA [9] , actually recognizing a specific target site on 
the 22-kD cx-zein gene [10]. The mutant gene greatly reduces 
storage zein protein, thus changing the endosperm texture. The 
02 gene has pleiotropic effects, and its function is complex. 
Proteomic analysis led Damerval to conclude that the 02 gene 
plays a key role in many metabolic pathways [11]. Microarray 
analysis indicated that 58 genes were up-regulated and 66 genes 
were down-regulated in o2 mutants [12]. Transcription profiling 
using GeneC ailing™ indicated that 70 genes corresponding to 23 
functional groups were up-regulated and 81 gene fragments 
belonging to 16 groups were down-regulated in o2 [13]. Some 
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studies reported differential expression of the o2 mutant gene in 
plants with different genetic backgrounds, resulting in varying 
levels of expressed oc-zeins [14] and lysine content [15,16]. In 
addition, the OPAQUE2 protein can interact with many other 
genes or proteins, such as the cyPPDK gene [17] and the ribosome- 
inactivating protein [18]. 

Pioneer Corporation [19] and the University of Missouri [20] 
exploited three SSR markers (phi 112, umcl066, and phi057) 
inside the 02 gene. The three markers can be used to select among 
the three genotypes, 0202, O2o2, and o2o2, in maize inbred lines 
constructed using molecular-assisted selection (MAS). 

Opaque-2 mutants have many poor agronomic characteristics, 
and QPM was developed to overcome these defects. To increase 
the lysine content of elite normal maize inbred lines, QPM inbred 
lines were used as donors to introgress the o2 mutant gene into 
elite Chinese normal maize inbred lines using SSR marker phi057. 
Nearly all of the selected lines presented opaque kernels and 
elevated lysine concentration. In our previous work, we found that 
some o2 near-isogenic lines (NILs) have a wild-type kernel 
phenotype with unchanged lysine content, such as Liao2345o2 
and Dan598o2 [21]. Surprisingly, another o2 NIL derived from 
liao2345 has an entirely different kernel phenotype with opaque 
endosperm that is identical to a typical o2 mutant individuals. In 
this paper, Liao2345o2 is referred to as Iiao2345/o2-l, and the 
second o2 NIL is referred to as liao2345/o.?-2. 

The 02 gene is pleiotropic, its function is complex, and some 
unknown mechanisms of this gene still need to be studied. We 
were interested in the reason why the two o2 NILs carry the same 
CDS but have completely different kernel phenotype. The object 
of this study was to reveal the mechanisms underlying the 
reversion of liao2345/ o2- 1 . 

Materials and Methods 

Plant materials 

Maize inbred line CA339 was used as a non-recurrent parent 
(donor), and maize inbred line liao2345 was used as a recurrent 
parent (receptor). CA339 is a QPM inbred line derived from 
pool33, which is from CIMMYT, and selected by Maize Research 
Center of Institute of Crop Sciences of CAAS; and liao2345 is a 
normal elite Chinese inbred line, cultivated by Maize Research 
Institute of Liaoning Academy of Agricultural Sciences. Two o2 
NILs, Iiao2345/o2-l and liao2345/o.2-2, were selected from the 
two parents by MAS [21]. The pedigree and detailed information 
refer to Table 1. In each generation, co-dominant SSR marker 
phi057 was used to select heterozygous (O2o2) individuals. The 
selected heterozygous individuals were self-pollinated and the 
recessive homozygous (o2o2) individuals were selected for further 
study. SSR marker phi057 was inside the 02 gene, according to 
the SSR marker phi057 the genotype of the two isolines was 
designated as o2o2. All of the materials were grown in Beijing 
during the summer of 2012. All plants were self-pollinated by 
hand, and the date and time were recorded. The endosperms were 
picked at approximately noon 1 8 days after pollination (DAP) and 
immediately frozen in liquid nitrogen. 

Kernel phenotype and total zein analysis 

Seeds were harvested after maturity, and dried in a greenhouse. 
The kernel appearance and characteristics were recorded using a 
camera (SONY cx700, P mode). Zeins were extracted according to 
the method described by Wu YR, with slight modification [22 
25]. For each sample, mature kernels from three ears were mixed 
evenly, and soaked in distilled water for 6 hours. The embryos and 
seed coats were removed and the endosperm was finely ground in 



liquid nitrogen using a mortar. Fifty milligrams of the resultant 
powder was precisely weighed and transferred to a 2-mL 
Eppendorf tube and 400 uL of extraction solution that contained 
70% ethanol/2% 2-mercaptoethanol (v/v)/28% double distilled 
water was added to the tube. The mixture was vortexed for 30 s 
and maintained on a shaking table for at least 6 hours. The 
mixture was then centrifuged at 13000 rpm (Eppendorf, Centri- 
fuge 5424R) for 15 min and 100 jiL of the supernatant was 
transferred to a new 1.5-mL tube. Next, 10 jiL of 10% SDS was 
added to the supernatant liquid and the liquid was dried in a 
drying oven at 50°C. One hundred microliters of distilled water 
was then added in order to resuspend the zeins. Finally, 8-jiL 
samples (equal to 1 mg endosperm powder) were loaded into 
polyacrylamide gel lanes, separated by a 4% stacking gel and 15% 
separation gel, and stained by Coomassie brilliant blue R250. 

Chromosome fragment introgression analysis in bin 7.01 

In order to analyze the introgression integrity of the o2 mutant 
gene, 18 SSR markers were selected in bin 7.01 on the 7th 
chromosome. All of the markers were selected from the database 
MaizeGDB (http://www.maizegdb.org/). 

DNA samples were extracted from leaves using the improved 
CTAB method [26] . The quality of the DNA was ascertained by 
0.8% agarose gel. The DNA concentration was measured by 
NANODROP 2000, then diluted to 100 ng/uL. The PCR 
reaction system contained 200 ng of genomic DNA, 1 x Taq 
Reaction Buffer, 0.2 raM of each dNTP, 0.25 mM of each primer, 
and 0.025 U/jiL of each Taq Polymerase, for a total volume of 
20 jiL. The touch-down reaction procedure was performed on a 
CI 000™ Thermal Cycler (BIO -RAD). The PCR cycling 
conditions were as follows: initial denaturation temperature of 
94°C for 5 min, followed by 10 cycles of 94°C for 30 s, primer 
annealing at 65°C for 30 s (every cycle decreased 1°C), and 
extension at 72°C for 45 s; this was followed by 25 cycles of 94°C 
for 30 s, primer annealing at 55°C for 30 s, and extension at 72°C 
for 45 s; and a final extension at 72°C for 10 min. After the PCR 
procedure, samples were maintained at 4°C until analysis. 
Electrophoresis and silver staining were performed as described 
in the CIMMYT laboratory manual protocols [27]. The results 
were analyzed using the software program GGT [28]. Bands 
consistent with the recurrent parent origin were recorded as "A," 
bands consistent with the non-recurrent parent origin were 
recorded as "B," and sites not marked by the markers were 
recorded as "U." 

Validation of promoter differences 

Several specific primer pairs were designed and used to 
determine the differences between two o2 NILs in the promoter 
region of the o2 gene. The primers were as follows: oligoFl: 5'- 
AAAGGTAGTAGAAC AC AC AC GGGC-3 ' , oligoF2 : 5 ' - 
C GTTGGTAAGGC ATGGTGTCTAGC A- 3 ' , oligoF3 : 5 ' - 
TGC AC GTAGTTC GTC GTC C AC AT- 3 ' , and oligo R: 5'- 
CGGTGGTAGCAGCTCCCAGAAGG-3'. All of the primers 
were designed to amplify the region that included the TATA box 
of the gene [7]. High-fidelity polymerase (KOD FX DNA 
Polymerase) was used to amplify the region, and the PCR reaction 
procedure was as follows: initial denaturation temperature was 
94°C for 5 min, followed by 35 cycles of 98°C for 10 s, primer 
annealing at 55 °C for 40 s, extension at 68°C for 80 s; and a final 
extension at 68°C for 10 min. The PCR products were maintained 
at 4°C until analysis. The products were separated on a 2.0% 
agarose gel in TAE buffer. 
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Table 1. Name and pedigree of maize lines. 





Name 


Pedigree 


Genotype 


Endosperm phenotype 


Lysine content% 


liao2345/o2-1 


((liao2345xCA339)BC 5 )F 7 


o2o2 


Hard and vitreous 


0.32[21] 


liao2345/o2-2 


((liao2345xCA339)BC 5 )F 7 


o2o2 


Floury and opaque 


0.49 


liao2345 


Tie7922xShen5003 


0202 


Hard and vitreous 


0.28[21] 


CA339 


pool33 


o2o2 


Semi-hard and semi-vitreous 


0.45 [44] 



doi:1 0.1 371/journal.pone.00851 59.t001 



Sequence analysis of the promoter and CDS region of the 
02 gene 

Three pairs of primers that were based on the 02 gene sequence 
published in the MaizeSequence database (http://www. 
maizesequence.org/index.html, GRMZM2G0 15534) were de- 
signed using Primer Premier 5 software. The primers are listed 
in Table 2 and shown in Figure 1 . High-fidelity polymerase (KOD 
FX DNA Polymerase) was used to amplify the 02 gene sequence. 
The cycling conditions were as follows: initial denaturation at 
94°G for 5 min, followed by 35 cycles of 98°C for 10 s, primer 
annealing at 55°C for 45 s, and extension at 68°C for 120 s with a 
final extension at 68°C for 10 min. The PCR products were 
maintained at 4°C until analysis. The products were separated on 
1.5% agarose gel in TAE buffer. The target genes were then 
purified and ligated into pEASY-J$\unt Cloning Vector (Beijing 
TransGen Biotech Co., Ltd.) using Escherichia coli strain DH5(X as a 
host. Plasmids were extracted and sequenced using the commonly 
used M13F and M13R sequencing primers. The sequences were 
blasted with the published sequence to ensure that the correct gene 
was ligated and all the sequences were aligned using DNAMAN 
version 6 software. And then the sequences were submitted to 
GenBank, the accession numbers provided by GenBank were 
KF831423, KF831424, KF831425, and KF831426. 

Analysis of 02 transcript abundance 

The following two pairs of specific primers that spanned the 
third and the fifth exon of 02 gene were designed for analysis of 
transcript abundance: 02RT1-F: 5 ' -TC AGGAATAATC - 
C AGTGC AGAA-3 ' , 02RT1-R: 5'-TCGACGT- 

TAGCGTCGTTGTA-3', 02RT2-F: 5 '-TAGCAGCATC AG- 
GAATAATC CAGTG-3', and 02RT2-R: 5'- 
CTTAGCTCTTAGGGTCTCCATGTCC-3 ' . The predicted 
PCR products were 358 bp and 409 bp. GAPDH gene was 
selected as the reference gene. The primers used for GAPDH were 
as follows: GAPDH-F: 5'-CCCTTCATCACCACGGACTAC-3', 
and GAPDH-R: 5'-AACCTTCTTGGCACCACCCT-3'. Endo- 
sperm collected 18 DAP was ground into fine powder and total 



RNA was extracted using TRIZol® agent (Invitrogen). An oligo 
d(T) primer kit was used to reverse transcribe the mRNA 
(Promega Corporation). A total of 2 Jig RNA was used to 
construct the first strain of cDNA. The reaction system was as 
follows: 4 uL of 25 mM MgCl 2 , 2 uL of RT 10 x Buffer, 2 uL of 
10 mM dNTP, 0.5 uL of 2500 U RNAE inhibitor, 0.7 uL of 
1500 U AMV RT, and 0.5 |ig of oligo d(T) primer, for a total 
volume of 20 jiL. The reaction cycle was 10 min at room 
temperature, 15 min at 42 °C, 5 min at 95 °C, and 5 min at 2°C. 
The reaction mixture was diluted with double distilled water to 
50% of its original concentration for PCR detection. 

Results 

Kernel characteristic and total zein analysis 

The two o2 NILs had different kernel phenotypes under normal 
and transmitted light (Figure 2). The kernels of liao2345/ o2-\ were 
hard and vitreous, similar to those of its recurrent parent liao2345, 
while the kernels of liao2345/o2-2 were soft and opaque. Kernel 
zein concentration and accumulation patterns, as determined by 
SDS-PAGE, are shown in Figure 3. The relative quantification of 
zein indicated that the zein concentration of liao2345/oJ?-2 was 
largely inhibited, while liao2345/oJ?-l had a zein concentration 
similar to that of liao2345 and larger than that of CA339. It was 
clear that the zein accumulation pattern was different among the 
four lines. The mutant o2 fragment of CA339 was introgressed into 
two o2 NILs. However, the kernel texture and zein accumulation 
pattern of these lines were completely different from one another. 
The line Iiao2345/o2-l had the same pattern as liao2345, whereas 
liao2345/oJ?-2 had a decreased zein concentration and a different 
pattern to that of liao2345 and CA339. Compared with liao2345, 
the o2 NIL liao2345/o2-2 had no 22-kDa cx-zein and a very low 
concentration of 19-kDa (X-zein. Both liao2345/o2-2 and CA339 
had lower concentrations of (X-zeins than liao2345 and the cx-zeins 
concentrations of CA339 were higher than that of liao2345/<x?-2. 
However CA339 possessed a higher 27-kDa y-zein concentration 
than the other three lines. 



Table 2. Primers for target gene fragments. 




Names 


Amplified Length 


Ta Opt °C 


Sequences 


SeqIF 


1373 bp 


55 


5' ATG G AG C ACGTC ATCTC A ATGG 3' 


SeqIR 






5' TAGTAGCAAGCCTCACGCCAA 3' 


Seq2F 


1227 bp 


57 


5' TAGCAGCATCAGGAATAATCCAGTG 3' 


Seq2R 






5' TCAGCGACGCCTGCAAATAA 3' 


Seq3F 


1419 bp 


57.1 


5' TCCTCTTTCCTTCTCTCAGACAGCG 3' 


Seq3R 






5' CGGTGGTAGCAGCTCCCAGAAGG 3' 


doi:1 0.1 371/journal.pone.00851 59.t002 
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SeqlF 



1 

lOObp 



A 

Seq3F 



A 

Seq3R 



SeqlR 
f 

■ Ml 



Seq2F 



A 

Seq2R 



4761 



Figure 1. Structure of the 02 gene and specific site of each sequence primer. 

doi:1 0.1 371 /journal.pone.00851 59.g001 



Validation of the o2 mutant gene between the two o2 
NILs and chromosome fragment analysis 

The SSR marker phi057 indicated that the two o2 NILs 
possessed the o2 mutant gene; however, they displayed completely 
different kernel phenotypes. We analyzed the chromosome 
segments in bin7.01 (Figure 4). Eighteen SSR markers were 
selected for examination of the two o2 NILs. Of these markers, five 
of them had no polymorphism or could not mark the 
corresponding genotype; five indicated that both NILs had 
segments from the recurrent parent (Figure 4, red region); five 
indicated that both o2 NILs had segments from the non-recurrent 
parent (Figure 4, blue region); and three markers, umc2160, 
phil 12, and umcl409, displayed polymorphisms between the two 
o2 NILs. The line liao2345/ o2-\ had fragments from the recurrent 
parent, liao2345, at the site of umc2160 and phil 12, while 
liao2345/o2-2 had fragments from the non-recurrent parent, 
GA339, at these two sites. At the site of umcl409, Iiao2345/o2-l 



had fragments from GA339 and liao2345/oi?-2 had fragments 
from liao2345. 

Three markers at position 122.4 were used: phi057 was a co- 
dominant marker that was located in the sixth exon, and was used 
to select the o2 NILs in this study; umcl066 was also a co- 
dominant marker, located in the first exon, but was not widely 
used for MAS because no polymorphisms were found among 
many materials; and phil 12 was located near the promoter. The 
results demonstrated that the o2 gene of Iiao2345/o2-l and 
liao2 345/o 2-2 were the same between the first exon and the sixth 
exon (Figure 5, umcl066 and phi057). However, Iiao2345/o2-l 
and liao2345/o2-2 displayed differences at the promoter region 
(Figure 5, phil 12): liao2345/ o2-l had the promoter of its recurrent 
parent liao2345, whereas liao2345/o2-2 had the promoter of its 
non-recurrent parent, GA339. We concluded that the two o2 NILs 
have different promoters. 
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Figure 2. Maize kernel phenotypes of liao2345/o2-1, Nao23455/a?-2, and Mao2345; o2o2, o2o2, 0202. A: Photographs of intact kernels 
taken under normal light on light box. B: Photographs of intact kernels taken with transmitted light. C: Photographs of decapped kernels taken under 
normal light. 

doi:1 0.1 371 /journal.pone.00851 59.g002 
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Figure 3. Maize kernel zein accumulation of liao2345/o2-1, 
Mao2345/o2-2, Nao2345, and CA339 in mature seeds, detected 
by SDS/PAGE. Total zein loaded in each lane is equal to 1 mg of 
mature seed endosperm, and the relative molecular weights of proteins 
are marked at the right side, the unit is kDa. 
doi:1 0.1 371 /journal.pone.00851 59.g003 

Validation of promoter differences 

To verify that Iiao2345/o2-l and liao2345/o2-2 had differences 
in the promoter region, several primers were used to examine the 
two o2 NILs. oligo Fl, oligo F2, oligo F3, and oligo R were 
designed to form three pairs of primers that were used to amplify 
the promoter region of liao2345/oi?-l and liao2345/o<2-2. The 
target DNA fragments of Iiao2345/o2-l were amplified, while 
those of liao2345/o<2-2 were not (Figure 6), confirming that 
Iiao2345/o2-l and liao2345/o,?-2 contained different promoters. 

Sequence analysis 

The primer pairs Seq3F and Seq3R also did not amplify the 
promoter regions of liao2345/ o2-2 and CA339. However, they 
could be used to amplify the promoter regions of Iiao2345/o2-l 
and liao2345. Moreover, the products were sequenced and aligned 
with the B73 sequence and the identity was 100%. Asymmetric 
interlaced PCR was used to determine the 5' flanking sequence of 
the liao2 345/ o2-2 and CA339 promoters. The mutant sequence 
was blasted in NCBI and had high homology with the Bg 
transposable element. It is the rbg transposable element that is 
inserted in the locus [38-41]. The 5' and 3' sequence of rbg and 
the 02 flanking sequence of liao2345/oJ?-2 and GA339 were 
cloned (Figure 7), and we verified that the rbg transposable element 
was inserted into the promoter regions of the o2 gene of liao2345/ 
o2-2 and CA339. And the partial sequence of rbg was submitted to 
GenBank, the accession number is KF831427. The insertion 
created an 8-bp direct repeat (GGCACAGC), and rbg had a 5-bp 
inverted repeat. The insertion was located between the TATA box 
and initiator codon, and was exactly 250 bp from the initiator 
codon (ATG). 
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Figure 4. Analysis of chromosome segments introgression of 
02 based on SSR markers mapping to bin7.01. The positions of 
SSR markers are based on a map from IBM2 2008 Neighbors 7. The red 
fragments (A) are of recurrent parent origin, blue fragments (B) are of 
donor origin, and gray fragments (U) represent markers that have no 
polymorphism between the two parents, and have no information for 
that genotype. 

doi:1 0.1 371/journal.pone.00851 59.g004 



The sequences of the 02 alleles, from the initiator codon (ATG) 
to the stop codon (TAG), were cloned and sequenced, including 
six exons and five introns. Figure 8 displayed multiple alignment of 
the 1 17 bp of the sixth exon. This alignment indicated that x 15544 
(supplied by GenBank) has -1 frameshift compare to the other six 
sequences. However, B73, Dhuang212, and liao2345 are wild-type 
inbred lines. Moreover, the B73 sequences were supplied by the 
database MaizeSequence (GRMZM2 GO 15534) [29]; therefore the 
stop codon should be a TAG, not a TGA. 
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Figure 5. Test of three SSR markers inside the 02 gene by PAGE, phil 12, umc1066, and phi057 can distinguish the three genotype of the 02 
gene. The amplified region of phil 1 2 located in the promoter of 02 the gene; and umd 066 in the first exon of the 02 gene; and phi057 in the sixth 
exon of the 02 gene. 
doi:1 0.1 371 /journal. pone.00851 59.g005 



Using the published sequence of B73, a wild-type inbred line, 
we deduced the OPAQUE2 protein sequence. Multiple alignment 
analysis (Figure 9) indicated that both of the o2 NILs acquired the 
donor's (CA339) coding sequence and were different from their 
recurrent parent (liao2345). In comparison to liao2345, liao2345/ 
o2-\ and liao2345/oJ?-2 had five amino acid deletions, six amino 
acids insertions, and 10 amino acid substitutions. Interestingly, 
there were only two sites at which amino acids substitutions 
occurred between the two o2 NILs. At these sites, a valine (V) was 
substituted by an alanine (A) and a glycine (G) was substituted by a 
valine (V). All of these amino acids are non-polar and hydropho- 
bic; consequently, they may not affect the structure of the 
OPAQUE2 protein. Further, in the OPAQUE2 protein, the 
leucine zipper domain of the two o2 NILs includes a substitution of 
a non-polar hydrophobic alanine (A) to a polar hydrophilic serine 
(S). The identity between B73 and liao2345 is 99.77%, whereas 
the identity between CA339 and liao2345 is 96.27%. Clearly, 
there are differences between the mutant gene and the wild-type 
gene, but the differences are not large and do not affect the 
function of OPAQUE2 protein. 
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Figure 6. Analysis of promoter difference between two o2 NILs 
by normal PCR on agarose gel. The size of the amplified region of 
three pairs primers are 582 bp, 648 bp, and 723 bp. This region contains 
the TATA box and initiation codon of 02 gene. 
doi:10.1371/journal.pone.0085159.g006 



Analysis of 02 transcript abundance 

Two pairs of specific primers were constructed to examine the 
transcript abundance of the 02 gene: both of them spanned the 
third and the fifth exon of the 02 gene. The transcript quantity of 
the 02 gene was obviously different between the two o2 NILs: 
liao2345, Iiao2345/o2-l had the same level of transcript, whereas 
the transcript of liao2345/o2-2 and CA339 were largely inhibited 
(Figure 10). The results demonstrated that the insertion of the rbg 
transposable element inhibited the transcription of the 02 gene. 
The transposon knocked down the expression of the 02 gene, and 
this was sufficient to induce a mutant phenotype. 

Discussion 

Detection of SSR marker phi057 indicated that two o2 NILs 
had successfully obtained the o2 mutant gene of CA339, so the two 
o2 NILs were designated as o2 mutant lines, and their genotype 
was o2o2. The germplasm of QPMs are associated with a recessive 
o2 mutant gene and modifier genes (o2 modifiers) [30,31]. These 
genes are responsible for the production of high lysine content and 
hard vitreous or semi-vitreous endosperm. The o2 modifiers can 
increase y-zein content by 2- to 3-fold [32,33] and the mechanism 
of this increase has been studied using RNAi [22,34,35]. In our 
study, the kernel phenotype (Figure 2) and the zein accumulation 
pattern (Figure 3) indicated that the 19- and 22-kDa ot-zeins were 
inhibited in liao2345/o2-2 and the 22-kDa (X-zein was almost 
undetectable, resulting in kernels that were opaque and soft. We 
speculated that the background of liao2345 contained no o2 
modifiers. Therefore, the maize inbred line, liao2345/ o2-2, is not a 
QPM inbred line. The maize inbred line Iiao2345/o2-l is also not 
a QPM inbred line because its lysine and y-zein concentrations 
were lower than those of GA339. In our previous study, liao2345/ 
o2-\ was determined to be o2 non-high lysine maize [21]. 

Chromosome fragment analysis demonstrated that the two o2 
NILs had different exchange patterns in bin 7.01. SSR markers 
umc2160 and phil 12 indicated that Iiao2345/o2-l had a large 
chromosome fragment of liao2345 (Figure 4). Recombination 
occurs during introgressing and backcrossing and progeny will 
have different recombination patterns. A very low frequency of 
recombination occurred inside the 02 gene, which means that the 
o2 gene of Iiao2345/o2-l was a fused gene, with a promoter 
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Inverted Repeats 
Direct Repeats 



agaagtttctgcagggaggcacagcCAGGGAAAACTTTAT ATAGAGTAAACCCTGggcacagcaagagagagagcccagca 

I II II 1 



02 



rbg 



02 



Figure 7. The 5' flanking sequence of the a? gene of Mao2345/o2-2. The capital letters represent the sequences of the rbg transposon, and 
the lower case letters represent the flanking sequences of the 02 gene. Red arrows indicate two inverted repeats of the rbg transposable element, 
and the blue arrows indicate two direct repeats, one of which was created by the insertion of rbg. 
doi:1 0.1 371 /journal.pone.00851 59.g007 



derived from liao2345 and CDS derived from CA339. However, 
liao2345/o2-2 inherited the whole mutant o2 gene from CA339. 
The results indicate that this may be the reason that two o2 NILs 
have different kernel phenotypes. 

In contrast to the published sequence (GenBank Accession No. 
XI 5544), our previous study indicated that the OPAQUE2 
protein of CA339 consisted of only 442 amino acid residues, 
shorter than that of XI 5544, due to a frameshift, resulting in the 
truncation of the coding sequence [21]. However, the newly 
published sequence of the B73 inbred line shows that there is no 
frame-shift mutation in CA339. In three wild-type inbred lines, 
B73, Dhuang212, and liao2345, the real stop codon is TAG 
instead of TGA (Figure 8). We retrieved another 47 partial 
sequences of the maize o2 gene from GenBank. Only two of these 
sequences were the same as x 15544 and the three partial 
sequences had -1 frameshift relative to the other sequences. The 
shifted partial sequences are NM 001111951.1, x 15544, and 
AY1 09364. We found that two o2 NILs, Iiao2345/o2-l and 
liao2345/oJ?-2, contained the same coding sequence of CA339. 
Between these two o2 NILs, there were two amino acid 
substitutions, valine (V) to alanine (A) and glycine (G) to valine 
(V). In both cases, the amino acids were non-polar and 
hydrophobic. Moreover, these amino acids are not in the basic 
motif of the OPAQUE2 protein and may not affect the function of 
the OPAQUE2 protein. Further study is needed to identify 
whether these changes disrupt the function of the OPAQUE2 



protein. A single amino acid substitution mutation from a non- 
polar hydrophobic alanine (A) to a polar hydrophilic serine (S) was 
found in the leucine zipper domain. This mutation might cause the 
OPAQUE2 protein to lose its spatial conformation, and thus, lose 
its function of transcriptional activation. However, there was no 
evidence as to whether this mutation might cause a loss of function 
in the OPAQUE2 protein. 

Transposable elements were first reported by McClintock [36— 
38] . Many transposable element systems have been characterized 
in maize. Transposable elements have many genetic effects, 
including mutation, chromosomal deletions, and rearrangements. 
The Bg-rbg transposable elements system is responsible for the 
instability of the mutable allele o2m(r) at the 02 locus [39,40]. Bg is 
an autonomous element and rbg is a non-autonomous element that 
is specific to the 02 gene. Nine o2 alleles were detected in the 
presence of an active autonomous Bg transposable element and 
five of these responded to the presence of Bg [41]. This 
phenomenon indicates that some o2 mutants were formed by the 
insertion of rbg transposable element. This type of mutation can be 
classified into a null mutation. In our study, we determined that 
liao2345/o2-2 and GA339 contained the rbg transposable element 
(Figure 7). The rbg has been cloned and sequenced [42]. It contains 
4500 nucleotide base pairs, making it easy to identify. Such a large 
fragment inserted in the promoter of the 02 gene could affect the 
transcript of the 02 gene, and we suspected that this insertion 
greatly inhibits the transcript abundance of the 02 gene. Our 



x 1 5544 CACCGCCGCCGCAG-ATTATGAGCTGCTGGGTCCAAATGGGGCCATACACATGGACATGT 

B73 CACCGCCGCCGCAGGATTATGAGCTGCTGGGTCCAAATGGGGCCATACACATGGACATGT 

Dhuang212 CACCGCCGCCGCAGGATTATGAGCTGCTGGGTCCAAATGGGGCCATACACATGGACATGT 

I i ao2345 CACCGCCGCCGCAGGATTATGAGCTGCTGGGTCCAAATGGGGCCATACACATGGACATGT 

CA339 CACCGCCGCCGCAGGATTATGAGCTGCTGGGTCCAAATGGGGCCATACACATGGACATGT 

I i ao2345/o2-1 CACCGCCGCCGCAGGATTATGAGCTGCTGGGTCCAAATGGGGCCATACACATGGACATGT 

I i ao2345/o2-2 CACCGCCGCCGCAGGATTATGAGCTGCTGGGTCCAAATGGGGCCATACACATGGACATGT 

x15544 ATTAGGCACTGCGGGTTTCGTGATCGCTGGGAACATTTTATTTGCAGGCGTCGCTGA 

B73 ATTAGGCACTGCGGGTTTCGTGATCGCTGGGAACATTTTATTTGCAGGCGTCGCTGA 

Dhuang212 ATTAGGCACTGCGGGTTTCGTGATCGCTGGGAACATTTTATTTGCAGGCGTCGCTGA 

I i ao2345 ATTAGGCACTGCGGGTTTCGTGATCGCTGGGAACATTTTATTTGCAGGCGTCGCTGA 

CA339 ATTAGGCACTGCGGGTTTCGTGATCGCTGGGAACATTTTATTTGCAGGCGTCGCTGA 

I i ao2345/o2-1 ATTAGGCACTGCGGGTTTCGTGATCGCTGGGAACATTTTATTTGCAGGCGTCGCTGA 

I i ao2345/o2-2 ATTAGGCACTGCGGGTTTCGTGATCGCTGGGAACATTTTATTTGCAGGCGTCGCTGA 



Figure 8. Multiple alignment of partial sequences of the sixth exon of the 02 gene. Sequences of xl 5544 and B73 are supplied by GenBank 
and MaizeSequence, respectively. The other five sequences are sequenced by us. Dhuang21 2 and liao2345 are wild-type inbred lines in China. Dashes 
indicate gaps and asterisks indicate identical residues. The stop codons are indicated by blue characters, and the insertions and deletions are 
indicated by red characters. 
doi:1 0.1 371 /journal.pone.00851 59.g008 
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B73 MEHV I SMEE I LGPFWELLPPPAPEPEPEREQPPVTG I VVGSV I DVAAAGHGHGGGDMMDQ 60 

I i ao2345 MEHV I SMEE I LGPFWELLPPPAPEPEPEREQPPVTG I VVGSV I DVAAAGHGHGGGDMMDQ 60 

CA339 MEHV I SMEE I LGPFWELLPPPAPEPE— REQPPVTG I VVGSV I DVAAAGHGDG— DMMDQ 56 

I i ao2345/o2-1 MEHV I SMEE I LGPFWELLPPPAPEPE— REQPPVTG I VVGSV I DVAAAGHGDG— DMMDQ 56 

I i ao2345/o2-2 MEHV I SMEE I LGPFWELLPPPAPEPE— REQPPVTG I VVGSV I DVAAAGHGDG— DMMDQ 56 



B73 QHATEWTFERLLEEEALTTSTPPPVVVVPNSCCSGALNVDRPP-VMEEAVMMAPAAVSSA 1 1 9 

I i ao2345 QHATEWTFERLLEEEALTTSTPPPVVVVPNSCCSGALNVDRPP-VMEEAVMMAPAAVSSA 1 1 9 

CA339 QHATEWTFERLLEEEALTTSTPPPAVVVPNSCCSGALNADRPPVVMEEAVTMAPAAVCSA 1 1 6 

I i ao2345/o2-1 QHATEWTFERLLEEEALTTSTPPPAVVVPNSCCSGALNADRPPVVMEEAVTMAPAAVCSA 1 1 6 

I i ao2345/o2-2 QHATEWTFERLLEEEALTTSTPPPAAVVPNSCCSGALNADRPPVVMEEAVTMAPAAVCSA 1 1 6 



B73 VVGDPMEYNA I LRRKLEEDLEAFKMWRAASSVVTSDQRSQGSNNHTGGSS I RNNPVQNKL 1 79 

I i ao2345 VVGDPMEYNA I LRRKLEEDLEAFKMWRAASSVVTSDQRSQGSNNHTGGSS I RNNPVQNKL 1 79 

CA339 VVGDPMEYNA I LRRKLEEDLEAFKMWRAASSVVTSDQRSQGSNNHTGGSS I RNNPVQNKL 1 76 

I i ao2345/o2-1 VVGDPMEYNA I LRRKLEEDLEAFKMWRAASSVVTSDQRSQGSNNHTGGSS I RNNPVQNKL 1 76 

I i ao2345/o2-2 VVGDPMEYNA I LRRKLEEDLEAFKMWRAASSVVTSDQRSQGSNNHTGGSS I RNNPVQNKL 1 76 



B73 


MNGEDP I NNNHAQTAG-LGVRLATSSSSRDPSPSDEDMDGEVE I LGFKMPTEI 


;rvrkrke 


238 


I i ao2345 


MNGEDP I NNNHAQTAG-LGVRLATSSSSRDPSPSDEDMDGEVE I LGFKMPTEI 


;rvrk- 


-ke 


237 


CA339 


MNGEDP I NNNHAQTAGGLVVRLATSSSSRDPSPSDEDMDGEVE I LGFKMPTEI 


;rvrk- 


-ke 


235 


I i ao2345/o2-1 


MNGEDP I NNNHAQTAGGI^RLATSSSSRDPSPSDEDMDGEVE I LGFKMPTEI 


;rvrk- 


-ke 


235 


I i ao2345/o2-2 


MNGEDP I NNNHAQTAGg[wRLATSSSSRDPSPSDEDMDGEVE I LGFKMPTEI 


;rvrk- 


-ke 


235 



B73 

I i ao2345 
CA339 

I iao2345/o2-1 
I i ao2345/o2-2 



SNRESARRSRYRK, 
SNRESARRSRYRK, 
SNRESARRSRYRK, 
SNRESARRSRYRK, 
SNRESARRSRYRK, 



AHLK LEDQVAQLKAENSCLLRR I AALI 
AHLK LEDQVAQLKAENSCLLRR I AALI 
AHLK LEDQVAQLKAENSCLLRR I ASLI 
AHLK LEDQVAQLKAENSCLLRR I ASLI 
AHLK LEDQVAQLKAENSCLLRR I ASLI 



3KYNDANVDNRVLRADME 298 

3KYNDANVDNRVLRADME 297 

3KYNDANVDNRVLRADME 295 

3KYNDANVDNRVLRADME 295 

3KYNDANVDNRVLRADME 295 



B73 TLRAKVKMGEDSLKRV I EMSSS-VPSSMP I SAPTPSSDAPVP PPP I RDS I VGYFSAT 354 

I i ao2345 TLRAKVKMGEDSLKRV I EMSSS-VPSSMP I SAPTPSSDAPVP PPP I RDS I VGYFSAT 353 

CA339 TLRAKVKMGEDSLKRV I EMSSSAVPSSMP I SAPTPSSDAPVPVPVPPP I RDS I VGYFSAP 355 

I i ao2345/o2-1 TLRAKVKMGEDSLKRV I EMSSSAVPSSMP I SAPTPSSDAPVPVPVPPP I RDS I VGYFSAP 355 

I i ao2345/o2-2 TLRAKVKMGEDSLKRV I EMSSSAVPSSMP I SAPTPSSDAPVPVPVPPP I RDS I VGYFSAP 355 



B73 AADDDASVGNGFLRLQAHQEPASMVVGGTLSATEMNRVAAATHCAGAMEL I QTAMGSMPP 41 4 

I i ao2345 AADDDASVGNGFLRLQAHQEPASMVVGGTLSATEMNRVAAATHCAGAMEL I QTAMGSMPP 41 3 

CA339 AADDDASVGNGFLRLQA-QEPASMVVGATLSATEMNRVAAATHCAGAMEL I HTAMGSMPP 41 4 

I iao2345/o2-1 AADDDASVGNGFLRLQA-QEPASMVVGATLSATEMNRVAAATHCAGAMEL I HTAMGSMPP 414 

I i ao2345/o2-2 AADDDASVGNGFLRLQA-QEPASMVVGATLSATEMNRVAAATHCAGAMEL I HTAMGSMPP 41 4 



B73 TSASGSTPPPQDYELLGPNGA I HMDMY 441 

I i ao2345 TSASGSTPPPQDYELLGPNGA I HMDMY 440 

CA339 TSASGSTPPPQDYELLGPNGA I HMDMY 441 

I iao2345/o2-1 TSASGSTPPPQDYELLGPNGA I HMDMY 441 

I i ao2345/o2-2 TSASGSTPPPQDYELLGPNGA I HMDMY 441 



Figure 9. Multiple alignment of deduced OPAQUE2 protein sequences. Dashes indicate gaps and asterisks indicate identical residues. The 
SSR site of umc1066 is underlined and the SSR site of phi057 is indicated by double underlining. The basic motif is marked with a red box, and the 
leucine zipper domain is marked by a yellow box. Blue ellipses indicate sites that differ between the two o2 NILs. 
doi:1 0.1 371 /journal. pone.00851 59.g009 



suspicion was then verified, as shown in Figure 10, the 02 gene of 
liao2345/oJ?-l and liao2345 were normally transcribed, while the 
transcript of liao2345/o2-2 and CA339 were inhibited to a large 
degree. We concluded that the two o2 NILs have the same coding 



sequence of GA339, and the two o2 NILs have different kernel 
phenotypes because they have different promoters: Iiao2345/o2-l 
has the wild-type promoter, while liao2345/o2-2 has the mutant 
promoter, as determined by their different patterns of recombi- 
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Figure 10. Examination of the 02 gene transcript abundance in developing endosperm (18 DAP) by RT-PCR. 02RT1 and 02RT2 are two 
specific primers spanning the third and the fifth exon of the 02 gene, the amplified region of ORT2 is larger than that of 02RT1; GAPDH is the 
reference gene. The templates of RT-PCR are labeled in each line, H 2 0 represents the blank control; RNA represents the negative control; cDNA is 
used to detect abundance of the 02 transcript. 
doi:1 0.1 371 /journal.pone.00851 59.g01 0 



nation. We deduced that the OPAQUE 2 protein of CA339 had 
the normal function of maintaining vitreous and hard endosperm. 
However, further studies, such as proteomic analyses, are required 
to determine whether the mutant OPAQUE2 protein has the 
same function as the wild-type OPAQUE2 protein. The unique 
maize inbred lines, Iiao2345/o2-l, liao2345/o2-2, and liao2345 
can help us to understand the role of the 02 gene, and its 
regulatory network in maize. 

We surmise that the reversion of liao2345/ o2-l to wild type was 
due to the recombination with the wild-type liao2345 promoter 
during introgression and backcrossing. Because the promoter 
sequences of liao2345 and Iiao2345/o2-l are identical. The 
promoter sequences may be different between the two lines, if it 
was due to the excision of the rbg element. The insertion of rbg in 
02 locus will create a 8-bp direct repeat, and the inserted base 
pairs may not be excised with the rbg element. Moreover, rbg is a 
nonautonomous element, and excision will not occur when the Bg 
element is not present in the background. The liao2345/o2-2 is 
very stable and did not revert; therefore, it contains no active Bg 
element. Only when reversion occurs prior to the first mitotic 
division of the primary triploid endosperm nucleus, will whole 
endosperm revertants (WER) be formed. Therefore, it is recom- 
bination that forms the reversion, not the excision of rbg. For the 
instability at the o2 locus [39], we need to use SSR marker phil 12 
when introgressing the o2 mutant gene into normal inbred lines to 
ensure that we have introgressed the mutant promoter. However, 
phil 12 is a dominant marker, and can not be used to distinguish 
homozygotes from heterozygotes. The results of the present study 



indicate that SSR marker phi057 is not a functional marker for the 
o2 gene. Even though there is polymorphism between receptors 
and donors in the simple sequence repeat region of phi057, the 
function of the 02 gene is not determined by this locus. Some wild- 
type maize lines, Duohuang29, Dan3130, and 9046, have no 
polymorphisms with CA339 [43], which indicates that they have 
the same sequence as CA339 at this locus. This further 
demonstrates that sequence variations at this site did not affect 
the function of OPAQUE2 protein. Our previous study demon- 
strated that there are two types of o2 mutant individuals: inbred 
lines CA339 and CA335 are different from Shandong2548 and 
Qi205, and have two different types of o2 mutant gene [44]. 
Functional markers should therefore be exploited according to 
specific o2 mutant genes, particularly SSR and CAPS markers, for 
MAS breeding of QPM. 
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